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Abstract 

^2 \ A concept of asymptotic symmetry is introduced which is based on a definition of 

symmetry as a reducibility property relative to a corresponding invariant ansatz. It is 
shown that the nonlocal Lorentz invariance of the free-particle Schrodinger equation, 
discovered by Fushchych and Segeda in 1977, can be extended to Galilei-invariant 
equations for free particles with arbitrary spin and, with our definition of asymptotic 
• symmetry, to many nonlinear Schrodinger equations. An important class of solutions 

of the free Schrodinger equation with improved smoothing properties is obtained. 

SI Introduction 

J> , It is well-known that the maximal Lie invariance algebra of the free linear Schrodinger 



X 



equation in three spatial dimensions 



Una i n 



m + 2^ A J ^ = °' (1) 

is the Schrodinger algebra sch(l,3), a Lie algebra which contains the Galilei and dilation 
algebras as well as some special conformal transformations [1-3]. It was quite surpris- 
ing therefore, when Fushchych and Segeda [4] showed that equation (1) is also invariant 
under an algebra of nonlocal pseudodifferential operators which is isomorphic to the Lie 
algebra so(l,3) of the three-dimensional homogeneous Lorentz group. 

The purpose of the present work is to establish similar results concerning nonlocal 
symmetry for some other equations of interest in mathematical physics. In particular, it 
will be shown that linear Schrodinger equations with linear and quadratic potentials with 
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arbitrary time-dependent coefficients and Hurley's Galilei-invariant wave equations [5], 
which describe free nonrelativistic quantum particles with arbitrary spin, are invariant 
under related algebras of nonlocal pseudodifferential operators. The results on nonlocal 
symmetries for linear Schrodinger equations with the potentials described above follows 
from the fact that these equations can be transformed to the free Schrodinger equation (1) 
and, therefore, inherit the nonlocal symmetry of the latter equation. The generators of the 
respective nonlocal Lie algebras for these Schrodinger equations necessarily have different 
representations than the generators of the algebras corresponding to equation (1), but the 
algebras associated, respectively, with these different equations are isomorphic. 

Analogous results to those for the linear Schrodinger equations mentioned above can 
be obtained for some nonlinear Schrodinger equations (NSEs). Specifically, a subclass of 
the family of NSEs proposed and discussed by Doebner and Goldin over the past sev- 
eral years [6, 7] in connection with representations of diffeomorphism groups and the 
corresponding Lie algebras of vector fields, as well as some related equations studied by 
Auberson and Sabatier [8] have been shown to be linearizable; i.e., they can be mapped 
by means of an appropriate change of variables to linear Schrodinger equations. In partic- 
ular, the potential in the latter equation can be chosen to be identically zero so that the 
solutions of these nonlinear equations can be related to the nonlocal invariant solutions of 
equation (1). These results make the situation for these linearizable NSEs analogous to the 
situation for linear Schrodinger equations with time-dependent linear and quadratic poten- 
tials, where the isomorphism of the respective Lie algebras was first proved by Niederer [9] 
for the linear Schrodinger equation with the usual (time-independent) harmonic oscilla- 
tor potential. In the cases examined by Niederer, however, the Lie algebras were those 
corresponding to the Schrodinger and oscillator groups whereas, in the present paper, we 
consider isomorphisms for nonlocal Lie algebras. 

Since the symmetries that we discuss are nonlocal, Lie's approach is not adequate to 
deal with them. A generalization of Lie's method, suitable for linear partial differential 
equations (PDEs), was suggested in [10] and is based upon the following commutator form 
of invar iance condition. 

Definition 1.1. An operator Q is a symmetry of a linear system of PDEs 



if and only if [L, Q]u = for each solution u of (2). 

Note that there is no restriction on the order of the operator Q (unlike the situation 
with Lie's methods [11]). In fact, Q need not be a differential operator, but may also be a 
pseudodifferential or integral operator. This non-Lie approach proved to be very effective, 
and wide classes of new symmetries were discovered for many equations of mathematical 
physics [12]. 

In the present paper we suggest an alternative symmetry criterion of invariance for 
nonlocal operators which is also suitable for nonlinear equations. 

Definition 1.2. We will say that an operator Q is a symmetry of a system of PDEs if 
and only if a corresponding ansatz, obtained as a solution of the equation 



Lu = 



(2) 



Qu = 0, 



(3) 



reduces the given system of PDEs. 
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Remark. When Q is a (linear) pseudo differential operator, one can construct an invariant 
ansatz (solution of (3)) by means of a Fourier transform, as shown in [13] (see also [3], 
Sec. 5.11). 

We prove that certain classes of nonlinear Schrodinger equations (NSEs) have a similar 
invariance when the time is larger than the squared modulus of the spatial variables in 
a certain well-defined sense. This "asymptotic symmetry" is based upon Definition 1.2 
when there does not exist an exact reduction of the given system of nonlinear PDEs by an 
invariant ansatz, but there does exist such a reduction in a well-defined asymptotic sense. 
(The precise definition will be given in Section 4.) This asymptotic symmetry situation 
is somewhat reminiscent of the situation encountered in scattering theory in which a 
solution of a nonfree equation approaches a solution of the corresponding free equation as 
t — > ±oo (cf. [14]). This scenario has been established for some NSEs with power-type 
nonlinearities [15]. In the present case, certain NSEs have a nonlocal asymptotic symmetry 
which is an exact symmetry of the free Schrodinger equation. 

As we shall discuss in more detail later, a number of authors have given Lie symmetry 
analyses for the NSEs that we consider [3, 11]. However, since the symmetries of interest in 
the present paper are nonlocal. Lie's algorithm cannot be applied to them and a different 
approach must be used to study them. 

The paper is organized in the following manner. In Section 2 we discuss the results 
for the free Schrodinger equation in further detail, extending the three-dimensional re- 
sults of [4] to any spatial dimension n > 2 and extending the nonlocal Lorentz-invariant 
solutions of the free Schrodinger equation to a larger class of symmetries. Nonlocal sym- 
metries of the type discussed in the present paper do not exist when n = 1. They are 
phenomena peculiar to higher spatial dimensions. Similar results are obtained for the 
three-dimensional Galilei-invariant free-particle wave equations of Hurley [5] for nonrela- 
tivistic free quantum particles with arbitrary spin. We have found a new representation 
of the Lorentz algebra for particles with arbitrary spin whose time evolution is described 
by these equations. In particular, the spin contributions to the angular momentum oper- 
ators do not satisfy the familiar angular momentum commutation relations, even though 
the total angular momentum operators satisfy the correct commutation relations with 
themselves and with the nonlocal generators. 

In Section 3 we show that, since linear Schrodinger equations with linear and quadratic 
(harmonic oscillator) potentials with arbitrary time-dependent coefficients can be trans- 
formed to the free Schrodinger equation (1), they inherit the nonlocal symmetry of the 
free equation. The corresponding nonlocal Lie algebras are isomorphic. Similar results 
are then discussed for the linearizable Doebner-Goldin and Auberson-Sabatier NSEs. 

In Section 4 we prove that several classes of NSEs possess nonlocal symmetries in an 
asymptotic sense when the time variable is sufficiently large. The collection of NSEs 
considered include several classes of standard type with power-type nonlinearities as well 
as the families of NSEs proposed and discussed by Doebner and Goldin [6, 7], Bialynicki- 
Birula and Mycielski [16], and Kostin [17]. Finally, Section 5 consists of some concluding 
remarks. 
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2 Free-particle equations 

We first discuss the exact nonlocal symmetry of the free-particle Schrodinger equation, and 
then consider analogous exact nonlocal symmetries for the Galilei-invariant wave equations 
of Hurley. 



Free-particle Schrodinger equations 

We define the standard operators po = ihdt, pj = —ihd X] (j = 1,2, ... ,n), and set 
PjPj 

Ls = Po — - where, in the latter definition, the convention of summing over repeated 
2m 

indices is used. Then the free-particle Schrodinger equation (1) takes the form Lsip = 0. 
Consider the operators (j, k = 1, 2, . . . , n > 2): 



and 



J jk = Xjp k - x k pj, (4a) 



Jo j = ^(f(p)G j + G j f{p)), (4b) 



where 

Gj = tpj — mxj, (4c) 

/„ \ 1/2 

P = I Yl PjPj ! an d / denotes a smooth function. 
\j=i J 

One easily verifies that the Schrodinger operator Ls is invariant under the operators (4) 
since 

[Ls,Jjk]=0=[L s ,Joj], j,k = 1,2, ... ,n. (5) 
The operators (4) satisfy the following commutation relations: 

[Jjhi Jqr] — ih(^5 r kJjq $qkJjr ^qjJkr ^jrJkq)i (6&) 

[Jjk, Jo q ] = ifr{8 qj J k - SkqJoj), (6b) 
[Joj,Jok] = -ih—Jjk, (6c) 

(j,k,q,r = 1,2,..., n > 2) where /' = f'(p) = — — . For f(p) = p, relations (6) are 

commutation relations of the Lotentz Lie algebra [4]. When f{p) ^ p, the invariance 
relations (5) are still valid. In this case, the operators {Jjk, Joq',j, k,q = 1,2, ... ,n} still 
form a Lie algebra isomorphic to the Lorentz algebra when / = \/p 2 + const. For other 
choices of /, the operators {Jjk, Jo q } do not form a Lie algebra although it is possible that 
they could be embedded in a Lie algebra of larger dimension than the Lorentz algebra. 
We will not consider this approach in the present paper. 
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Using the symmetries (5), one can derive fundamental solutions ip of Lgip = which 
are invariant under the operators (4) by requiring that Jjktp = 0, Joqip = (j, k,q = 
1,2,..., n). In terms of the Fourier transform ip of ip , this procedure leads to: 

V;(fc,t) = c/(fc)- 1 / 2 exp(-|^), (7) 

/„ y/ 2 

where k = ^ kjkj and c denotes a complex constant. This result was first obtained 

V=i / 

(for the case f(k) = k) in [13]. In order to calculate the inverse Fourier transform of (7), 
it is convenient to specify the function /. A convenient set of choices is: 

Case I: f(k) = k a , < a < 2n. 

This reduces to the Lorentz case when a = 1. One finds from (7) (see [18] for a 
derivation in the case a = 1): 



ip(x,r) = (2irh)~% / exp ( % -k ■ x ) ^{k,t)d n k 



(8) 

( m \2 4 1 V2 4/ 



\2mht) r(f) lFl \2 4'2' 2M 



n an im\x\ 



(valid in all dimensions n>2if0<a< 2n) in terms of the confluent hypergeometric 
function iF 1: and we have chosen the constant c in (7) so that (8) reduces to the standard 
Galilei-invariant fundamental solution 



/ m \ f / im\x l2 



when a = 0. For the case n = 3 and a = 1, one may use the recurrence relations for the 
confluent functions and the well-known relations 

1F1 (v + ^2v + l;2iz^j =T{l + u)e iz (^~ U J v (z) 
to write (8) in terms of Bessel functions (see [13, 3]): 



, , x i47T4/m\| / im\x\ 2 \ . .1 [ / m\x\ 2 \ . T /m|x| 2 

^ = {^ihi) exp v^rj N ' [ J - 4 l^rj +zj t ( 



4fa 



(10) 



Similarly, when a = n > 2 the same relation between iiq, and J;, given above can be used 
to write (8) in the following Bessel function form: 

it / m \- I im\x\ \ I mix \ \ V4 21 ( m\x\ 



ip(x, t) = -3— — — exp — — — - J 



' n 1 



\2nihtJ \ 4ht J V 8fit / 4-5 ^ 4fa 

where the following well-known relation between T functions has also been used: 



2 / 2 2z ~ 1 
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We note that, when n is a multiple of 4, the Bessel function in the above expression can 
be further simplified to a combination of sine, cosine, and algebraic functions. 

The fact that the fundamental solutions (8) have a slower time decay, t~(~2~4) than 
the corresponding Galilean result (9), which is seen to be t~~, was discussed in [18] for 
the case a = 1. The point was made that, although this property makes (8) unacceptable 
for use in quantum mechanics (when used as a kernel analogously to the usual use of the 
Galilean fundamental solution (9)) because the standard probability interpretation for 
the wave functions of free nonrelativistic particles is not obtained; its use may actually be 
more desirable than (9) for mathematical applications because of the smoothing properties 
which the linear operators that are constructed from them possess. (Analogous operators 
constructed from (9) do not have such smoothing properties.) 

One advantage of considering the cases a / 1 of Case I is that the corresponding 
convolution mappings formed from (8) or (10) have improved smoothing properties as 
a increases. Thus, generalizing the discussion in [18] for the case a = 1, we define the 
mappings G n : g -> G n (a)g: 

(G n (a)g)(x, t) = (27r/i)-S exp (^k ■ x"j (|£f )~ f exp (-^) flf(fc) dk, 

Qt 

and deduce that the maps G n (a) are smoothing in the sense that G n g have — (distribution) 

derivatives if g G L 2 (R n ) (n > 2). 

The above argument shows that the smoothing properties of mappings on L 2 con- 
structed with the fundamental solutions (8) or (10) increase as a increases. However, 
since f(p) = p a of Case I only increases algebraically with a, and a is bounded above 
by 2n, the smoothing properties are limited. This suggests that one consider functions / 
of exponential type in the variable k 2 : 

Case II: f(k) = exp (2(3\k\ 2 ) , (n > 2). 

In this case, the corresponding mappings G n have much improved smoothing properties 
because G n : L 2 (W 1 ) -> L q (W 1 ) for all q G [2,oo). In addition, it turns out that the 
corresponding fundamental solutions (i.e., the inverse Fourier transforms of (7)) have the 
same asymptotic time decay as t — > +oo as the Galilean fundamental solutions (9): 

( it\~% ( Ixl 2 

i/>(x,t)=c (2h(3+-) exp - ' 

The fact that the functions (11) have the same asymptotic decay as the Galilean funda- 
mental solutions (9) is explained by the fact that the former correspond to imaginary time 
translations of the latter: ip(x, t) = ip g (x, t — 2mih(5) with the choices c = (2tt?i)~z for the 
constant in (11). The increase in smoothness of the mappings G n associated with (11) 
relative to those associated with (9) is intimately connected with the fact that the time 
translations leading from (9) to (11) are imaginary. 

Free-particle wave equations for arbitrary spin 

According to the principles of quantum mechanics, a wave function which represents a 

1 3 

free nonrelativistic particle with spin s = — , 1, — , . . . should have 2s + 1 components. The 
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equations obtained by Hurley [5]: 

/ ihdths+i 
1 



L H ?p = 



ns hs 



\ 



2mhs 
i 

2mhs 



S • V hs+l 02s+l,2s-l 

K-p 2s -l,2s+l hs-1 J 



x 
n 



o. 



(12) 



are a system of 6s + 1 equations in which the first row in (12) gives the equations of motion 
for the (2s + l)-component wave function ip and the remaining two rows are equations 
of constraint defining the redundant functions x (2s + 1 components) and (2s — 1 
components) in terms of ip. Sj and Kj (j = 1,2,3) are spin matrices with dimensions 
(2s + 1) x (2s — 1), (2s — 1) x (2s + 1), respectively, and can be chosen to satisfy the 
following relations [5] (with K* the adjoint of K and Eijk the completely antisymmetric 
Levi-Civita symbol): 



SiSj + K*Kj 



ishEijkSk + h 2 s 2 5ij. 



(13) 



The symbol I m denotes the m-dimensional unit matrix, and a 5 denotes the zero matrix 
with a rows and b columns. Equations (12) reduce to the equations derived by Levy- 

Leblond [19] when s = - and by Hagen [20] when s = 1. They are Galilei invariant 
by construction. For convenience in the demonstration of nonlocal invariance of equa- 
tions (12) in the discussion to follow, we have used a different normalization than Hurley. 

By substituting the equations of constraint (the second and third rows of (12)) into 
the equations of motion for ip (first row of (12)), and using (13), one finds that each 
component of if) satisfies equation (1). In order to establish the nonlocal symmetry of 
equations (12), we define new representations of the commutation relations (6). Thus, in 
place of the operators (4), we define (with j, k = 1,2, 3): 



Jjk = ( x jPk - XkPj)ks+i - -^( x jPk - hPj), 



and 



J 0j = —(f(p)G J +G J f(p)), 



where 



Gj = (tpj - mxj)I 6s+ i + Xj, 



(14a) 



(14b) 



(14c) 



with 



/ 02s+l,2s+l 02 S +I,2s+1 02 S +I,2s-1 \ 
02s+l,2s+l 02 S +I,2s-1 
02s-l,2s+l 02s-l,2s-l 



-s, 

2s 3 
2s 3 



(14d) 



Then, using (13), one verifies that solutions of (12) are invariant under the operators (14) 
in the sense that [L#, $]ip = when Lfftp = 0, with d = Jjk, Joj (j, k = 1,2, 3) and that 
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{Jjk, Joj] h k = 1,2,3} satisfy the commutation relations (6) with {Jjk, Joj} replaced by 
{Jjk,Joj}, respectively. 

We can now construct solutions of (12) in the form 

4 = column (^-^S-P^-^-K-p^ (15) 

with ip(x,t) = P(t)F(x,t), where (3(t) has 2s + 1 components and F(x,t) is a (scalar) 
solution of (1). Then, inserting (15) into (12) and using (13), we infer that the components 
of (3 must be constant. In particular, F(x,t) can be taken as the invariant solution (10). 
However, in order to write the solutions (15) in complete detail, one needs representations 
for the spin matrices Sj and Kj (j = 1, 2, 3). Examples of these are given in [5]. 

Finally, we note that the operators Joj defined in (14b)-(14d) are pseudodifferential 

operators and that the terms {^jPk — ^kPj) (j,k = 1,2,3) in (14a), which play the 

role of spin angular momenta in the present case, are noncanonical in the sense that they 
do not satisfy the usual angular momentum commutation relations. Thus, they do not 
generate the usual representations of the three-dimensional rotation group corresponding 
to spin s. Nevertheless, as we noted above, the operators {Jjk, Joj', 3, k = 1,2,3} satisfy 
the correct commutation relations. 



3 Nonfree equations possessing exact nonlocal Lorentz sym- 
metries 

In this section we discuss how the exact nonlocal symmetry of free-particle equations, 
discussed in the preceding section, can be extended to some nonfree equations which 
describe interactions between particles. The existence of nonlocal symmetries for these 
equations follows from the fact that appropriate Lie algebras can be constructed which 

are isomorphic to the Lie algebra formed by the operators in equations (6) with = 1. 

P 

We discuss this situation for two cases: (1) linear Schrodinger equations with linear and 
quadratic potentials (with arbitrary time-dependent coefficients), and (2) some classes of 
nonlinear Schrodinger equations. 

Linear Schrodinger equations with linear and quadratic potentials 

The symmetry properties of these equations have been well-studied, especially in one space 
dimension. We refer to a recent discussion of the latter case in the context of coherent 
states and squeezed states [21], from which many references may be traced. 

Our approach to these equations is based on the fact that they can be transformed 
to the free Schrodinger equation (1). The results are analogous to those obtained by 
Niederer [9] for linear Schrodinger equations with time-independent harmonic oscillator 
potentials which showed that the Lie algebras of the oscillator and Schrodinger groups 
are isomorphic. Our transformation results generalize the treatment of Niederer and are 
a direct extension of those of Truax [22], Bluman [23], and of Bluman and Shtelen [24] 
for the case of one spatial dimension. However, our objective is different than that of the 
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authors cited above in that we are interested in algebras of nonlocal symmetries rather 
than in algebras of point transformations. 
For Schrodinger equations of the form: 



2 



ihd t ip(x,t) + — Ai/;(x,t) = (a(t)\x\ 2 + bj{t) Xj +c(t))4>(x,t), (16) 

n 

with \x\ 2 = x j x j for n > 1, where a, bj (j = 1,2, ... ,n), and c are arbitrary functions 

3=1 

of t, we define the following transformation: 

4>(x,t) = exp (A{t)\x\ 2 + Bj{t)xj + C{t))^j u(y,r), (17a) 

with 

T W = / ^V)^, yj = a{t)x j + p j (t), (17b) 

where a{t), A{t), C(t), Pj(t), and Bj(t) (j = 1,2, ... ,n) are to be expressed in terms of 
a(t), bj(t) (j = 1,2, . . . , n), and c(t). By substitution of (17) into (16), one can choose the 
former coefficients so that (16) reduces to equation (1) for the function u{y,r): 



2m 

1 



ihd T u(y,T) + —A y u(y,T) = 0. 



For the standard harmonic oscillator: a = ^moJ 2 , bj = (j = 1, 2, . . . , n), c = 0; (17) 
reduces to the transformation obtained by Niederer: 



ip(x,t) = (sec(ut)) 2 exp y — tan(cjt)|x| j u(y,r), (18) 

with yj = sec(ut) X j (j = 1,2, ... ,n) and r = — tan(cji). By taking u(y,r) to be the 

u> 

Galilean fundamental solution (9) expressed in the y, r variables, one obtains the standard 
fundamental solution for the harmonic oscillator ([9] and [25], p. 63). In a similar manner, 
we can use the mapping (18) to obtain nonlocal invariant solutions by taking for u(y,r) 
the nonlocal invariant solution (8) of equation (1): 

(, n a -p I n a\ / • I |2 \ 

m n 1 U ~ 4 ) v n a n im \y\ \ ^ 
2^Tr) ^|T lFl U-4 ; 2 ; ^rJ' 

and obtain 

t fmojtan(iot)\-^ T (f - § ) 



muj 



2mhsm{ujt)) \ 2irih j T (§) 



imoj j , \ i 1 2 " l n a n im ^\x 



x exp -—— tan(wt)|x| \Fi 



2h v n 1 J \2 4' 2' hsm(2ut) 



2 
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Remark. An explicit form for the generators of the Lorentz algebra for the harmonic 
oscillator is as follows: 



Joj = (pGj + Gjp) , P = Y^PiPi , n>2, 



where 

pj = cos(u>t)pj — miosin(tot)xj, 

and 

Gj = — sin(cot)pj — mu> cos(ut)xj, j = 1, 2, . . . , n, 

in terms of the same coodinates Xj, and momentum operators pj used in Section 2. We 
use the same angular momentum operators (j, k = 1,2, ... ,n) as defined in (4a) and 
obtain the following commutation relations in place of (6b), (6c) for the free case: 

[Jjk, Joq] = ifr{5qjJok ~ dqkJoj), 

[Joj,Jok] = -ihJjk, j, k = 1,2, ...,n > 2. 

We derive similar results for linear potentials by setting a = 0, c = 0, bj / (j = 
1, 2, . . . , n > 1) and obtain mappings of solutions of (1) onto solutions of (16): 

/ i ( |fr| 2 t 3 \\ 

ip(x,t) = exp ( -- ItbjXj + J 

at 2 

yj = aXj + —bj, (j = 1,2, ... ,n) T = a 2 t. 
a = const; and mappings of nonlocal invariant fundamental solutions (with n > 2): 
m \f-f r(f-f) (n an im\x\ 2 it i\b\H z 

u(y ' T) ~ VhMH) r(§) 11 V 2 " I ; 2 ; ~2hT + 2h bjXj + ~8^n~ 

nX,t) ~ r (f ) P V & V 5) 6m J J 

n an im|x| 2 it i|6| 2 t 3 



XiFi ---; -; rr~ + ^rbjXj + 



2 4 2 ' 2ht 2h JJ 8mh 

Similar results can be given for the case of a harmonic oscillator driven by an external 

force fj(t) (j = 1, 2, . . . , n), for which a = -mu 2 = const and bj = —fj, whose Galilean 

fundamental solutions were given in [25]. However, since the expressions obtained are 
somewhat unwieldy, we shall not give those results here. 
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Nonlinear Schrodinger equations with exact nonlocal symmetry 

Doebner and Goldin [6, 7] considered a class of nonlinear Schrodinger equations which 
were suggested by their studies of dissipative quantum theory based on group theoretic 
considerations relating to groups of diffeomorphisms on Euclidean spaces and the corre- 
sponding Lie algebras (algebras of vector fields). Related equations were considered earlier 
by Sabatier [26] and subsequently by Auberson and Sabatier [8] and by Auberson [27]. 

We first consider a subclass of the Doebner-Goldin (DG) equations which are lin- 
earizable in the sense that they can be mapped to linear Schrodinger equations by point 
transformations. The nonlocal symmetry of the linearizable DG equations follows from 
the isomorphy of the Lie algebras of these equations with the nonlocal Lie algebras of 
linear Schrodinger equations. 

The DG equations can be written in the form: 

ihik = ( - J- A + V(x, t)U + —Raty, V0V> + KD 1 £ CjRjty, (19) 
\ m / j=1 

where D and D' denote constant diffusion coefficients, and the real-valued nonlinear 
functionals Rj(ip,ip) (j = 1,2,3,4,5) are given by: 



(20) 



p p p 

where p = tpip, and j is related to the usual probability current density j by 

j = — j = ^(V'VV' - ipVip). 

Lie symmetry analyses of these equations have been discussed by several authors [28-30]. 
References to discussions of analogous equations by other authors can be found in [7]. 

It was pointed out by Doebner and Goldin that the subfamily of their equations (19), 
(20) defined by the following relations between the coefficients D, D', Cj (j = 1, 2, 3, 4, 5): 

D = D' Cl = -D'c 4 ; D'(c 2 + 2c 5 ) = D'c 3 = 0, (21) 

are linearizable and that the corresponding solutions can be constructed from solutions of 
linear Schrodinger equations by means of "nonlinear gauge transformations" (see also [28] ) . 
The nonlinear gauge transformations of DG are given by: 

^ ft = N(ip) = |V>|exp^[7ln|V>| + AArg^, (22) 

with 7, A real numbers (and A / 0). Doebner and Goldin show that, given the rela- 
tions (21) among the coefficients, if ip is a solution of (19), (20); then tp' = N(ip) is a 
solution of the following linear Schrodinger equation: 



ih ,_ _h 2 
~A ~ ~2^A 2 



ti = -z=i*W + V{x,tW, (23) 
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when 



4m . 4m 2 Z) 5 
A = 1 —D 'c 2 



and 



2mDA 

7 = -— , (24) 



■ , , , 4m . Am 2 D 2 , _,. . . . , 

provided that — — D C2 H 775 — < 1. Since the gauge transformations satisfy the group 

n n 2 

law 7V 71 Ai • N r2 a 2 = ^ 7l+ Ai72 AiA 2 > t ne gauge transformations inverse to (22) are given 
by: 

if, = N~ l {ij') = exp [-A-^ln^'l + A _1 Arg^]) . (25) 

This mapping is analogous to the mapping (17) which transforms solutions of the free 
Schrodinger equation to solutions of the linear Schrodinger equation with linear or quad- 
ratic potentials. 

For the case when the potential V is identically zero, ip' is a solution of the free 
Schrodinger equation with the mass m replaced by the "effective mass" raA. For ex- 
ample, we may take the solution analogous to (8): 

^ {x ^-[^nht) r(f) lFl V2"4'2'^^J (26) 

for < a < 2n and n > 2. More generally, one can also consider mappings from solutions 
of the DG equations (19), (20) to solutions of a linear Schrodinger equation with one 
of the potentials discussed in the first part of this section and then map to solutions of 
the free Schrodinger equation by using the results discussed there. The composition of 
this sequence of mappings yields transformations of solutions of (19), (20) to solutions of 
free Schrodinger equations without the necessity of assuming that the linear potential in 
the DG equations is identically zero. 

Auberson and Sabatier [8] (AS) considered the following NSE (for convenience, we set 
h = 1 and 2m = 1): 

# t (x, t) = (-A + V)1>(x, t) + s^i>(x, t), (27) 

where s is a real parameter. For s < 1 AS use the following linearization transformation: 
ip = \ip\ex.p(-i9), t = (l-s)-h', 0(x,t) = (1- s)h'(x,t'), (28) 

which transforms equation (27) to the following linear Schrodinger equation: 
iip' t ,(x,t') = -Aip'(x,t') + (1 - s)~ 1 y(x)V' / (x,t') 

for the quantity 

ip'{x,t') = ip (x,(l - syh'^j exp(-i0 / (x,t'))- (29) 
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AS also linearize (27) when s > 1, but the linear equations thereby obtained are not 
Schrodinger equations so we shall not discuss them. 

If we consider invariant solutions analogous to (26) for (29): 



then the corresponding invariant solutions of (27) are obtained from (28). 

As a consequence of the analogies between the inverse gauge transformation (25) of 
the DG equations or of the linearization transformation (29) of the AS equation with the 
transformation (17) between solutions of the free Schrodinger equation and solutions of 
linear Schrodinger equations with linear or quadratic potentials, we see that the symmetry 
algebras corresponding, respectively, to the free Schrodinger and to the linearizable DG 
or AS equations are isomorphic. 

4 Asymptotic symmetry 

For general nonlinear Schrodinger equations, one does not expect the nonlocal symmetries 
for the equations discussed in the preceding sections to extend in an exact form because 
those results depended on the fact that the equations were either linear or linearizable. 
Lie's algorithm [3, 11] gives a general method for the investigation of local symmetries 
and solutions of differential equations, including nonlinear ones; and this method has 
been applied to many types of equations, both linear and nonlinear. However, as we 
have noted in the Introduction, since the symmetries that we are discussing are nonlocal 
and are defined in terms of pseudodifferential, rather than differential, operators; Lie's 
approach and related techniques are not adequate to deal with them. Because of this 
difficulty of extending nonlocal symmetries to general nonlinear equations, we propose to 
use a definition of symmetry based on a reducibility property (see Definition 1.2 and the 
discussion below). When there is no exact reducibility, we introduce a weaker concept of 
asymptotic symmetry. 

Power-type nonlinearities 

Consider an NSE of the form 



= 1,2, ... , n), where the nonlinear function F depends in general on the solution ip, its 
complex conjugate ip, derivatives of these functions through the second order and, unless 
a statement is made to the contrary, we set h = 1 in the present section. 

Following Definition 1.2, we will say that an operator Q is a symmetry of equation (30) 
if and only if the corresponding ansatz, obtained as a solution of (3), reduces (30) to 
a system of PDEs in fewer independent variables or, as a limiting case, to a system 
of ordinary differential equations (ODEs). This approach is especially useful when one 
wants to extend symmetries of a system of linear equations (in the present case the free 
Schrodinger equation (1)) to a system of nonlinear equations. 





(30) 
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By extending an argument in [13] for Case I with a = 1, we conclude that the following 
ansatz is invariant under the algebra (4): 

ij(x,t)=<j ) (t)g(t,x), (31) 

where (j) is an arbitrary function of t and g has the form (8) (or (10) in the special 
case when n = 3 and a = 1). According to Definition 1.2, we say that equation (30) is 
invariant under the algebra (4) if and only if the ansatz (31) reduces (30) to an ODE for the 
function (f)(t). In the following definition, we introduce a concept of asymptotic symmetry 
when an exact reduction does not exist, but a reduction does exist in an asymptotic sense. 

Definition 4.1. We will say that equation (30) has the asymptotic symmetry (4) if 
and only if the ansatz (31) reduces (30) to an ODE for <f>(t) in the asymptotic region 
m\x\ 2 <C 2t. 

In this section we will first discuss NSEs with power-type nonlinearities and then con- 
sider several cases of derivative nonlinearities. We first treat the case in which the non- 
linear term in (30) is of the form: 

F = A(V^) fc V, (32) 

where k denotes a positive real number (not necessarily an integer) and A denotes a 

complex (coupling) constant. Lie symmetry analyses of equations of the form (30), (32) 

have been discussed by many authors (see [3] for a summary). In addition, many authors 

have investigated the existence of solutions to such equations in various Banach and Hilbert 

In 

spaces. For dimensions n > 3, many of these results require that < k < because 

n — 2 

the proofs use the Sobolev embedding theorem. See [15] for a summary. Our results are 
not subject to this restriction. 

To investigate the asymptotic symmetry of (30), (32), we look for a solution of the 
form (31) where g(x,t) is a solution of the free Schrodinger equation defined by (8): 



I -A ( 771 \ 2"4 I U ~~ l) Ei I 71 a U lm \ X \ \ 

9{x > f) = (2^) 4%r 11 V 2 " 4 ; 2 5 —) 



(33) 



with < a < 2n for spatial dimensions n > 2. (If n = 3 and a = 1 we may, of course, 
use (10) instead.) Then, substituting these expressions into (30) and (32), assuming that 
4> depends only on t, and using the fact that g satisfies the free Schrodinger equation 

idtg = Ao, we obtain the following equation for (f>: 

2m 

<j>t = -i\(M k (gg) k ^, k > 0. (34) 

Since <\> is assumed to depend only on t, the above derivation is only consistent if the 
quantities (gg) k are independent of the spatial coordinates xj (j = 1,2, ... , n). This is 
not true in general, so we consider the limit of large t or, more precisely, values of the 
variables Xj (j = 1, 2, . . . , n) and t such that m\x\ 2 <C 2t. Then, using the small-argument 
expansions 

1 F 1 (a; C; ,) = l + -, + -L- T i- + ... (35) 
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for the confluent functions, we obtain the following asymptotic result for gg: 
99 



SHraf 1 + U^' 1 '' <36> 




and (34) becomes, to leading order in the quantity 



m\x\ 2 



2t 

(\ 2k 
r(f-f) 
r(f) 



= -*M^7 (37) 



From the form of equation (37), we find that its solution must have the form <p(t) = 
(3 exp(—ih(t)) with (5 a real constant and h(t) a real-valued function of t. Substitution of 
this expression into (37) yields an equation for h(t) which has the following solution: (uJ 
a real constant) 

/ \ 2k 

M(HA(£) (»-^r(|g)^ (38) 

where Xk(t) = t~( n ~ z) k+l when k ^ {n — and = Int when k = — . The 

asymptotic solution to (30), (32) is now obtained by substitution into (31): 

ip(x,t) = 0exp(-ih(t))g(x,t) 

with g given by (33) and h given by (38). 

We note that the asymptotic result for |^| = (ipip)? corresponds (apart from the con- 
stant P) to the asymptotic result for the solution g(x, t) of the free Schrodinger equation, 

1 (ip\ 

whereas the asymptotic value of Arg = — In I = j contains effects of the nonlinear terms 
(32). 

Similar results can also be obtained when the nonlinear term in (30) is a linear combi- 
nation of power-type nonlinearities such as, for example: 

F = —ao^ — aifyip)^ — a,2(ipip) 2 ip, (39) 

where aj, (j = 0, 1,2) are real constants, and 02 / 0. Lie symmetries of equations (30), 
(39) were discussed by Gagnon and Winternitz [31]. 



The Doebner-Goldin and related equations 

We next consider asymptotic symmetry results for the DG equation (19), (20) and then 
discuss the relationship between these results and some others for these equations. 

To show that the DG equations have asymptotic symmetry in the sense of Definition 4.1, 
we follow the procedure used for NSEs with power-type nonlinities and look for a solution 
ip of (19), (20) of the form (31) with g(x,t) defined by (33) or (10) with the appropriate 
powers of h again inserted. We obtain the following linear equation for ip by virtue of the 
homogeneity property of the nonlinear functionals Rj (j = 1,2,3,4,5): 

' Pi D ^ 

ihfo = V i p+ i -^R 2 (g,9)ct> + nD'Y,CjR J (9,m- (40) 

3=1 
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Using the derivative relations — \Fi(a; c;z) = — iiq(a+l; c+1; z) and the small-argument 

dz c 

expansions (35) for the confluent functions, we obtain the following asymptotic results for 
the Rj(g,g) (j = 1, 2, 3, 4, 5): 




|2 



2ht 




^=( 1 -s)^(^) +0 

_ , _. / a \ 2 2a 2m / m|x| 2 \ z / ^ / m|x| 2 \\ 



n 2 (ra + 2) 2 ht \ 2Ht J \ \ 2fit 

Thus, we see that the DG equations are asymptotically invariant in the sense of Defi- 
nition 4.1 if we consider the case of an identically zero potential V and omit the Rj 
functional with j = 2, 3, 4, 5. Then, solution of (40) gives the following asymptotic result 
(m\x\ 2 < 2ht): 

(/ (y. \ vrun \ 
-iD' Cl (1 - — J — ln(t) J g(x, t) (41) 

with k a complex constant. 

For reasons of consistency, we must show that a solution ip of (19), (20) obtained as 
in (25) from a solution ip' of (23) is consistent with our asymptotic symmetry result (41). 
This can be done by noting that, for the case when the potential V is identically zero, ip' 
is a solution of the free Schrodinger equation with the mass m replaced by the "effective 
mass" mA as in (26). Using the small-argument expansion (35) and related expansions for 
the arctangent and logarithmic functions that occur in Arg ip' and ln(|^'|), respectively, 
we obtain from (25): 



, , , / mA \ § r (f - f ) / 2miD , (v (f - § ) / mA \ 




(42) 



itt (n a.\ ( a \ 

exp -xd-sJ+H 1 -^) 



m|x| 2 ' 
2/it 



We will compare this expression with the asymptotic symmetry result (41), which can be 
written in the form: 



n a 



( rah \ 2 4 r m - f) / r . / a \ mn , N 
t) <* « ( — j exp (-i (l - -) - In(t) + 



/n a 



)]) 



exp 



a \ m|x 



2n) 2ht 
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The expressions (42), (43) will be compared in the spacetime domain m|x| 2 <C 2ht by 
setting 



k = /3exp(— id) with and 5 real. 
One then obtains 







A 2 4 



(44) 



(45) 



and the coefficients of In t in the exponents of (42) and (43) agree because of the equality 
D = D'ci, which is part of the conditions (21) of DG required for linearizability. Equating 
the constant parts of the phases of (42) and (43) gives: 



/ a 
= n 1 — — 
V 2n 



In) 



1) 



mD 



In 



mA 
2^h 



2mD 



h 



In 



r(f) 



(46) 



Thus, the solutions of the DG equations (19), (20) obtained by their gauge equivalence to 
solutions of the linear Schrodinger equation (23) with V = is consistent with the asymp- 
totic symmetry result (41) provided that the constant k is chosen to satisfy (44)-(46). 



Equations of Kostin and of Bialynicki-Birula and Mycielski 

We consider NSEs of the form: 



(47) 



where and £2 are real constants. The first logarithmic term ln (l^l 2 ) ip was originally 
proposed by Bialynicki-Birula and Mycielski [16] and the second logarithmic term, which 

involves the phase of ip: (2i) _1 ln ( = ) , was first proposed by Kostin [17] in connection 

\ipj 

with studies of dissipation effects in quantum mechanics. It is appropriate to discuss equa- 
tion (47) in this section because Doebner and Goldin have shown [7] that their equation 
(19), (20) extends to include the nonlinear terms in (47) when the parameters 7 and A in 
the gauge transformation (22) are time-dependent. Symmetry analyses of equation (47) 
with £2 = have been investigated (cf. [ 32, 3]) and existence results for this equation 
(with £2 = 0) when an appropriate class of linear potentials is also included have been 
summarized by Cazenave [15]. 

Looking for solutions of (47) in the form (31), (33); we obtain the following equation 
for <j): 

iHt = 6 In (H 2 | 5 | 2 ) <j> - I ln (|0 4>- (48) 
Then, assuming that m|x| 2 <C 2ht and using the expansions (35), we obtain (36) and 



434 



W.W. Zachary and V.M. Shtelen 



Equation (48) becomes, to leading order in the quantity 



m x 



+ 



2h 



In = 



+ 



2ht ' 

2h 



(50) 



where the last term on the right-hand side comes from the logarithm in (48) and vanishes 
— — — is an even integer. From the form of eq.(50), we find that its solution must 

have the form (f>{t) = (3 exp(z<5(i)) with (3 a real constant and 5(t) a real-valued function. 
Substitution of this expression into (50) yields the following equation for 6(t): 



<5 t = -|ln[ I 



2f rn xn-§ / r(f-«4) \ 2 \ M 6/n_^ 

\27rfa/ ^ r(§) ) y h ^ hU 8 J 



(51) 



The solution of this equation has different forms according as £2 is zero or nonzero. 
Case A. £ 2 = 0. 

2 N 



In /3 2 



■p / n a 

1 I 2 4 

r(§) 




a 



+ ( n - - ) ( In 



/ m \ 



+ 1 



(52) 



x exp (^ (n — ^ tint + i(J , ( = real const. 

Tl OL 

Case B. £ 2 / and — — — / an even integer 

In this case, equation (51) can be written in the form: 

I h (-t ) *>) 4 (» - !) — (-¥) + ¥ (M) - (-¥ 




/^M"-f f r(f-f) 



r(f) 



ex P( -- 



Integrating this equation from i* > to +00, and using integration by parts for the lni 
term, we obtain 



(§) 



I (->«•-( 



n a 
4 ~ 8 



and 



(f)(t) = (3exp(i5(t)) = (3 exp ( i 

.a 



x exp 



6 



n — 



2"/ 



lni 



i 

.a 

s 6 



1 In {?(—T % f r(t " f) 

2 r v2vrJ rm 



n 



exp 



6* 



£1 



(§) 
ft 




(53) 
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roo 

in terms of the exponential integrals E\(x) = — Ei (— x) = I exp(— t)t~ l dt (where the 

J X 

integral is understood as a principal value integral when x < 0). In Case B, the terms 

— — — ) 7r are absent when — — - is an even integer. The final 
4 8/ 24 

I 1 2 

vn \x\ 

asymptotic solutions to (47) are obtained to leading order in the quantity — - — by substi- 
tuting (52) and (53) into the equation ip(x,t) = <f>(t)g(x,t), where g(x,t) is given by (33). 

5 Concluding remarks 

We have shown that nonlocal Lorentz symmetry, previously known to be valid for free- 
particle Schrodinger equations [4, 13], is also valid in a modified form for the Galilei- 
invariant linear wave equations of Hurley [5], which describe the time evolution of free 
quantum particles with arbitrary spin. In Section 3 we have discussed similar results for 
linear Schrodinger equations with linear and harmonic oscillator potentials with arbitrary 
time-dependent coefficients. For nonlinear equations, we have shown that a subset of the 
nonlinear Schrodinger equations introduced by Doebner and Goldin [6, 7] in connection 
with studies of dissipative effects in quantum mechanics as well as some of the related 
equations discussed by Auberson and Sabatier [26, 8] have exact forms of these symmetries. 
Moreover, we have also shown that several classes of nonlinear Schrodinger equations 
- those with power-type nonlinearities as well as the full set of equations proposed by 
Doebner and Goldin - have asymptotic nonlocal symmetry in the sense described in the 
present paper. 

The question of the nonlocal symmetry of nontrivial (i.e., nonfree) multiparticle Schro- 
dinger equations is open. In general, we expect that the concept of nonlocal symmetry 
may be helpful in the analysis of such equations. 
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